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We present a powerful and general approach to describe the coupling of Majorana fermions to 
external leads, of interacting or non-interacting electrons. Our picture has the Klein factors of 
bosonization appearing as extra Majoranas hybridizing with the physical ones. We demonstrate the 
power of this approach by solving a highly nontrivial SO(Af) Kondo problem arising in topological 
superconductors with M Majorana-lead couplings, allowing for arbitrary M and for conduction 
electron interactions. We find that these topological Kondo problems give rise to robust non-Fermi 
liquid behavior, even for Fermi liquid leads, and to a quantum phase transition between insulating 
and Kondo regimes when the leads form Luttinger liquids. In particular, for M = 4 we find a long 
sought-after stable realization of the two-channel Kondo fixed point. 

PACS numbers: 73.23.-b,74.78.Na,72.10.Fk,03.67.Lx 



One of the most influential recent discoveries in con- 
densed matter physics is that topological phases support- 
ing Majorana fermions can be engineered in heterostruc- 
tures based on s-wave superconductors and materials 
with strong spin-orbit coupling [l|. This breakthrough 
development, transforming so- far elusive ideas Q for Ma- 
jorana based fault-tolerant quantum computers into a po- 
tentially feasible perspective, drives immense theoretical 
and experimental activities (see Ref. for reviews). 

Majorana fermions are localized, robust, zero-energy 
excitations. Testing the zero bias anomalies^ Q arising 
from electrons tunneling onto them forms one of the main 
experimental directions with promising results so farQ. 
Transport can also inform on the Majoranas' quantum 
computational potential, if it indicates that the nonlo- 
cal "topological qubits", formed of pairs of them obey 
quantum dynamics 

In all transport problems, a key role is played by the 
coupling of Majorana fermions 7 to leads of conduction 
electrons ip. Here we introduce a general picture of this 
coupling (Fig. [J), observing that vb = T/ can be bro- 
ken into Majorana (r) and charge density (/) parts us- 
ing bosonization. Upon T hybridizing with the physical 
Majorana, the charge sector couplings get organized in 
a simple structure. Rendering the charge sector trans- 
parent is a key virtue of our picture. It allows one to 
approach qualitatively new problems involving the inter- 
play of Majorana fermions and electron interactions, be 
those between conduction electrons or due to the charg- 
ing of the superconductor. 

We demonstrate this by solving a highly nontrivial 
problem that includes both types of interactions. This is 
the topological Kondo effect[7|, which, for M Majorana- 
lead couplings implements a novel, SO(Af) Kondo prob- 
lem. In Ref.[Z|, we could solve this for the simplest, SO(3) 
~ SU(2) case for noninteracting leads. Our picture lets 
us explore vastly more general settings, allowing for ar- 
bitary M and interacting conduction electrons. We find 
a number of striking features. These include stable, non- 



FIG. 1: A Majorana mode 7 on a superconductor heterostruc- 
ture coupled to a lead of conduction electrons. The Majorana 
r arises from the conduction electrons ip = Tf. It hybridizes 
with 7 while / changes the charge density (shown with waves). 

Fermi liquid (NFL) behavior, even for Fermi liquid leads; 
a quantum phase transition between Kondo and insulat- 
ing regimes due to the competition between the Kondo 
effect and the suppression of electron tunneling in Lut- 
tinger liquids; and a long sought-after, stable realization 
of the two-channel Kondo fixed point for M = 4. Our pic- 
ture delivers these results with ease through mapping the 
problem to a quantum Brownian motion (QBM) model. 

We begin by setting up our Majorana-lead coupling 
picture more explicitly. We consider a superconducting 
structure with M to t localized Majoranas jj . Their quasi- 
particle operators obeyQ 



7j=7]i {7j!7fc} = 2<5. 



(1) 



Ordinary fermionic modes Cjk=^('jj+ijk) arise from pairs 
of Majoranas; in a fermionic system M tot is thus even. 

We will work with half infinite (x > 0), single- 
channel leads furnishing effectively spinless conduction 
electrons. This can be achieved in several Majorana 
realisations [lj, including the nanowire based setups of 
recent experiments [fij]. Bosonization [8| transforms the 
Hamiltonian of lead j into a quadratic problem, 

fx f 

Ho(<PjA) = / dxK{d x e 3 ? + K-\d xVj ) 2 , (2) 

even if interactions are present. Here 0j,<Pj are bosonic 
fields encoding the charge density pj = % <Pl . They obey 

[<p(x),6(y)}=4iriG(x-y), (3) 
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where 0(x) is the Heaviside function, and satisfy^ 
(p(0) = (d x 8)(0) — at the endpoint of the lead. Electron 
interactions enter Eq. ^ through the Luttinger parame- 
ter K and through rcnormalizing the velocity u. We have 
K < 1 (K > 1) in the repulsive (attractive) and K = 1 
in the noninteracting regimes. 

Working at energies much below the gap, M leads cou- 
ple to the superconductor through M Major anas: 
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>«/2 



+ H.c, 



(4) 



where tj is the tunneling amplitude, ipj annihilates elec- 
trons at the end of lead j, and the phase exponential 
e ±*x/2 changes the number of electrons on the supercon- 
ductor by ±1. 

Bosonization fractionalizes the electron operator as 
ipj(x) = Tjfj(x), whereQ fj = -^e 10 ^ 2 at the endpoint, 
up to a numerical factor to be absorbed in tj and with 
being a the short distance cutoff. The Klein factors Tj 
are often omitted but for us they are crucial. In ensuring 
that ipj ipki li anticommute, they extend Eq. ([!]): 



{r i! r fc } = 25 



Ok, 



(5) 



This allows one to view Tj as additional Majorana 
fermions which can, for example, form ordinary fermion 
modes dj = \{~fj + iTj) with the physical Majoranas. It 
also means that the terms in Eq. (Q| 
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M 
.7 = 1 



e i9j/2 e ix/2 



H.c. 
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with real ti,t[, S. Rotating from 71, 72, Ti to new Majo- 
ranas 71.2.3 decouples 73 to give 



H t = i^ 2 \/S 2 + ^ + !Lj£ cos 2 (0l/2 ). 
a a 



(8) 



This immediately recovers two key results of Ref. [|§ ex- 
panding the root for t\ w tf (which holds in the nontopo- 
logical case) we find that the leading perturbation, which 
is also the most relevant in the sense of the renormaliza- 
tion group (RG), is Cooper pair tunneling cx cos 2 (#i/2); 
and that for t\ ^> t' 2 , aS 2 near the topological regime the 
S and t'i terms give the same perturbation. 

We now move on to the topological Kondo effect 0). 
In this Kondo effect strong correlations emerge even for 
noninteracting leads from the interplay of the charging 
energy E c of the superconducting island (which is now 
connected to ground by a capacitor) hosting the Majo- 
ranas and the ground state degeneracy associated with 
themQ. This degeneracy can be understood in terms 
of the fermion modes Cjk' they have zero energy (up to 
exponentially small corrections in the Majoranas' sepa- 
ration) which leads to a 2 Mtot_1 fold degenerate ground 
state, where the (—1) in the exponent is because ^2 c ]k c jk 
mod 2 is tied to the electron number parity of the island. 

The fermion modes Cjk provide the topological qubits 
and the Kondo effect arises from coupling these to con- 
duction electrons via the Majoranas. For low tempera- 
tures, voltages and weak coupling, T, V, \tj\ <C E c , the 
physics is dominated by virtual transitions connecting 
the lowest energy charge state of the island to the neigh- 
boring ones with ±1 extra electrons. These are captured 
by the effective Hamiltonian 



factorize into Majorana-Majorana and charge sector 
parts. This is our main observation. 

Before applying this to the topological Kondo problem, 
we warm up with the setup of Ref. |5f tunneling between 
a superconductor and a Luttinger liquid lead (M =1). 
Recovering the couplings obtained there is a good valid- 
ity check of our picture, which also allows us to highlight 
some of its advantages. The superconductor is grounded 
in this setup, thus we can drop e 1 */ 2 for now. If the sys- 
tem is in the topological regime the coupling is through a 
single Majorana, which enters the problem only through 
7iTi. This is related to the parity djdi=i(l + i7iri) and 
acts as <7i in the two state space of a c\j fermion (in the 

(— l) CljClj =<73 basis). This immediately recovers the cou- 
pling of Ref. [a circumventing the need for the augmented 
Jordan- Wigncr transformation employed there. 

In the nontopological regime, there is another Majo- 
rana 72 which couples to the lead and 71. The phase of 
C12 and ip can always be chosen so that the coupling is 

^ cos(0i/2) , „ sin(0i/2) 



Hr.tr — 
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j^k=i 
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where \f k 



tjtk/E c 



> 0. Eq. 



is obtained by a 



Schricffcr-Wolf transformation keeping the leading order 
terms in tj/E c . As shown in Refs. [7J, the nontrivial lead- 
qubit couplings of the first term implement an SO (A/) 
Kondo problem (with a spinor "impurity spin" through 
7,-7fc). Ref. solved Eq. © for SO(3)~SU(2), applying 
the Affleck-Ludwig conformal field theory method[10(|. 
This is however unsuited for interacting leads and be- 
comes complicated for M > 3. As we now show, our 
picture handles both challenges with ease. 

Substituting ipj = Tjfj into Eq. (O, we again find that 
Majoranas enter through jjTj Importantly, jjTj 

commute with each other and thus can be diagonalizcd 
simultaneously, jjTj = ±i, where we can absorb the sign 
in 9j. We get 



-i0 k /2 iBj/2 



j 



(10) 
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Eq. pop is the central formula underlying all our sub- 
sequent analysis. We reduced the Kondo problem into 
one for the charge densities, eliminating the Majorana 
degrees of freedom through the 7j Tj hybridizations. 

For A - = 0, Eq. ([TU| is known to map to QBM which, 
for A+ fc = A + < 0, K > 1, is related to the M > 3 channel 
Kondo modeld [H [H. Our SO(M) problem is in the 
A + > sector, and we have K < 1 in its lead-Majorana 
implementation. Eq. ([10]), for M = 3, A+ fe = A+, A^ = 0, 
also appeared as an unphysical description of quantum 
wire tri-junctions, used for illustrating the dangers of 
omitting Klein factors 0- For us Eq. (|10l) arises precisely 
from using Klein factors correctly. 

In what follows, we apply the RG to our full prob- 
lem, focusing on K < 1 and allowing for the anisotropy 
(X^ k 7^ A + ) and \T. terms arising in a physical realiza- 
tion. The key aspects of the physics will be shown to 
come from the X^ k term, extending the relation to QBM 
to this more general case. The results will be used to 
obtain the low temperature behavior of the Kubo con- 
ductance Gki between leads k,l. The setup (for M = 5) 
is sketched in Fig. [2] . 

The weak X^ k flow under rescaling a — > ae l is obtained 
from the operator product expansion, 



A'o = vq • if [with vq 
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k ' 
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Here v is the density of states of the leads. The cou- 
plings Xjj do not renormalize. As the A^, terms transfer 
charge between the leads [see Eq. (|5J)], the first term in 
Eq. (jlip promotes a suppression, the second term an en- 
hancement of tunneling processes. This corresponds to 
the competition of two distinct mechanisms: the suppres- 
sion of electron tunneling in Luttingcr liquids jl4|. and 
the enhancement of the coupling due to the Kondo ef- 
fect. Depending on which of these wins one will find a 
markedly different low temperature behavior, described 
by the decoupled lead and the Kondo fixed points. 

The transition between the two cases is governed by a 
repulsive, isotropic (A^~ fc = A + ) intermediate fixed point 
at v\ + * = 2K(m-2) ' w i*h exponent — 1 in the relevant, 
isotropic direction. (The fixed point is attractive from 
the orthogonal directions.) This is consistent with the 
weak coupling regime for K < 1, but we believe that 
the existence of an isotropic fixed point governing the 
transition remains true in general. 

For A < A + *, where A is the typical bare value of A^., 
the low temperature behavior is dictated by the decou- 
pled fixed point. In the opposite case, Aj" fc flow to large 
values while becoming more and more isotropic (consis- 
tently with the flow around A + *). The crossover to strong 
coupling is at the scale of the Kondo temperature, Tk ~ 
E c e~ 1 / uX . The analysis of the strong coupling regime 
becomes transparent after rotating 9 = (#i,...,#m), 

= (ipi, . . . , (Pm) to decompose them to Rq = vq • 6, 



components rj, kj along the directions orthogonal to Vo- 
This canonically decouples the Kq,Rq sector, while the 
r,k sector has /^2Ho(kj,rj) perturbed by 



i(w fc -Wj)-r/2 ,\., 

H eS = -y a+ y -^w, • d x k, 



(12) 



where w,- • w/ 



jt, implying that the minima of 



the first term form a. (hyper)triangular lattice. Near the 
Kondo fixed point A + — > oo, r tends to be pinned to one 
of these minima. This, at the same time, suppresses^ 
9 x k, thus the A~ term can be neglected in this regime. 
The Xjj couplings, therefore, only give small, marginal 
perturbations (through the Rq, Kq sector) again, show- 
ing that the QBM model indeed captures the essential 
physics in both the weak and strong coupling regimes. 

Using QBM results [HI we immediately find that the 
leading perturbations at the Kondo fixed point, given 
by r tunneling between the adjacent minima, have di- 
mension A — 2K ^ [ ~ 1 ^ > ■ The Kondo fixed point is thus 



stable as long as K > 



> 



it is robust against 



2(A/-1) ' 2' 

weak conduction electron interactions and its robustness 
is enhanced as M increases. Due to their non-integer di- 
mension, these charge conserving processes do not admit 
a free fermion description. The SO(A/) Kondo problems 
thus give rise to NFL behavior even with Fermi liquid 
(K = 1) leads. Such stable NFL fixed points are absent 
from conventional Kondo systems 15[ . Our method al- 
lows us to prove that this remarkable feature is a generic 
property of the topological Kondo effect, vastly general- 
izing the M = 3, K = 1 result of Ref. 0. 

In particular, for M = 4, K = 1 we find A = 3/2, 
as for the two-channel Kondo fixed point [loj]. Through 
relabeling ipj one can indeed map our problem to the 
two-channel Kondo model as A^ fc — > A + , A^ — > 0. The 
topological Kondo effect thus provides a long sought-after 
stable realization of this fixed point that does not hinge 
on (but is robust a gain st) having NFL leads, in contrast 
to earlier proposals [16|. Placing the two channel Kondo 
in the QBM context, we also find a new theoretical per- 
spective, alternative to Refs. 13, 17. 18. 

We can now apply our findings to the Kubo conduc- 
tance Gki- The phase diagram in terms of A and K 
is sketched in Fig. [5J Its topology is dictated by the 
QBm[I3|, but the physical meaning of the phases is spe- 
cific to the topological Kondo effect. Tuning A or if, the 
system undergoes a quantum phase transition, switching 
Gk^Li from to at T — > 0. The transition using A 

is especially appealing, as A is gate tunable, in principle, 
in the nanowire realizations of recent experiments @. 

Near the decoupled lead fixed point, Gki vanishes as 
G u ~ T 2 / R - 2 as T — > 0. This is the known suppression 
of electron tunneling between Luttinger liquid leads fl4|. 
with the exponent coming from the first term of (|11[) . 
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FIG. 2: The phase diagram of the topological Kondo effect 
in terms of the typical bare A + coupling A and the Luttinger 
parameter K. The arrows indicate the flow of the couplings 
under the RG, and Gk^i is the low temperature Kubo con- 
ductance. The curved line separates insulating (A + — > 0) and 
Kondo (A + — > oo) phases. The inset shows the sketch of the 
setup for M = 5. The central rectangle is an s-wave su- 
perconductor and the light bars between the Majoranas are 
semiconductor nanowires. 

For T <C Tk near the Kondo fixed point, we have 
G k i = ^T-($ki - j I ) + a kl T 2A - 2 , where a M are nonuni- 
versal constants. The fixed point (T — > 0) value fol- 
lows from the emergent boundary conditions r = fixed, 
d x \i = at A + — > oo, and can be obtained from an im- 
mediate generalization of the calculations of Ref. 0. Note 
that Gu = x^' ^ stable Kondo fixed point (A > 1) 

2 

thus comes with Gu violating the %- limit for single- 
channel normal conduction. This is due to the emergent 
boundary conditions translating into correlated, multi- 
particle Andrcev processes by which holes, not only elec- 
trons can be backscattered[9|. The T 2A ~ 2 dependence is 
due to second order corrections in processes tunneling r 
between the minima in Eq. (fTiZj) . (The first order correc- 
tions of the current-current correlators underlying Gui 
vanish [HI].) The convergence to an enhanced conduc- 
tance through such nontrivial power laws gives a clear 
signature of the NFL Kondo physics. 

In conclusion, we have introduced a bosonization based 
picture for Majorana-lead couplings relevant to ongoing 
transport experiments. The key feature is the appear- 
ance of Klein factors Tj, virtually extending the num- 
ber of Majoranas in the system. We have illustrated 
the utility of our approach applying it to two-terminal 
normal-superconductor transport with interacting leads 
and solving the SO(M) topological Kondo problem. We 
expect that our picture will be a useful starting point 
for a number of new problems exploring the interplay of 
Majorana fermions and strong correlations. 
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